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Abstract: We check the validity of the widely used classical apsidal motion formula as 
a function of orbital parameters, stellar structure, and stellar rotation rate by comparing 
dynamical calculations of the periastron advance with the static tidal formula. We find 
that the classical formula gives very accurate results when the periods of the low order 
quadrupole g, f and p modes are smaller than the periastron passage time by a factor of 
about 7 or more. However, when this condition is not satisfied, the difference between the 
classical formula and the exact result can be quite large, and even periastron recession can 
result. The largest difference arises when one of the low order modes of the star is nearly 
resonant with an integer multiple of the orbital frequency minus twice the rotation rate 
of the star. The resonance of higher order g-modes (number of radial nodes > 4) with 
the orbit is very unlikely to cause significant deviation from the classical result because 
of their weak coupling to the tidal force and thus their small contribution to the apsidal 
motion. Resonances involving rotational modes of the star are also unlikely to make much 
contribution to the apsidal motion because of their small overlap with the tidal force, even 
though they have periods comparable to the periastron ge time. 

We apply our work to two famous binary systems (AS Cam and DI Her) which show 
abnormally small apsidal motion, and conclude that dynamical effects are unimportant for 
these systems, i.e. the static tide assumption is an excellent approximation. 
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1. Introduction 



When the gravitational field of a star diff"ers from that of a Newtonian point mass, 
the orbit of its companion will deviate from a Keplerian orbit. To lowest order, a per- 
turbation to the 1/r gravitational potential causes the periastron to rotate. This is the 
origin of apsidal motion. There are primarily three effects which cause deviation from 
the 1/r gravitational potential: the General Relativistic correction to Newtonian gravity, 
the quadrupole moment that arises due to the rotational distortion of a star, and the 
quadrupole moment due to tidal distortion. The first two effects are relatively easy to 
calculate and are well understood. The third effect, the modification of the gravitational 
potential due to tidal distortion, displays more complex behavior and is the subject of this 
paper. The derivation of the classical apsidal motion formula, due to tidal distortion, was 
first worked out by Cowling (1938) and Sterne (1939) and is widely used in the interpreta- 
tion of observed apsidal motion. Their derivation follows from the approximation that the 
star can instantaneously adjust its equilibrium shape to the time-dependent tidal force as 
it orbits around the center of mass of the system. This approximation is clearly valid in 
the limit that the periods of the low order modes of oscillation of the star are much smaller 
than the orbital period (for orbits of low eccentricity) or the periastron passage time (for 
eccentric orbits). However, this approximation is expected to fail when the mode periods 
are of order the relevant time scale. Indeed, as Papaloizou and Pringle (1978) have pointed 
out, resonances can significantly modify the classical formula. 

In view of the wide use of the classical apsidal formula, it is important to understand 
where the classical approximation is valid and where it fails. In this paper, we model the 
tide dynamically as a linear superposition of the non-radial oscillations of a star excited by 
the gravitational tidal force of the primary. We calculate, using a fourier series expansion 
of the oscillation amplitude, the quadrupole gravitational potential perturbation due to 
these oscillations, and the resulting periastron shift. Our result is compared with the 
classical apsidal motion formula for stars of different polytropic index and mass, and for 
various orbital parameters. The rotation of the star and its modification of resonances 
is also considered. We also consider the effects of rotational modes on apsidal motion. 
However the modification of p and g mode eigenfunctions due to stellar rotation, which 
can be quite substantial, is not included in our calculations. 

2 



The plan of our paper is as follows. In §2 we present the basic equations for the 
perturbation to the gravitational potential arising from dynamic tides and discuss the 
validity of the classical apsidal formula for nearly circular orbits. In §3, we examine the 
effects of tidal resonances on apsidal motion as a function of modal and orbital parameters. 
In §4, we compare our results for polytropic stars with the classical apsidal formula to find 
the system parameters for which deviation from the classical prediction is expected. The 
effect of rotational modes on apsidal motion is considered in §5. In §6 we examine the 
apsidal motion of two peculiar binary systems (AS Cam and DI Her) which have observed 
periastron advances that differ significantly from the predicted values. We investigate 
whether the discrepancy might in fact be due to the inapplicability of the classical apsidal 
formula. Our results are summarized and discussed in §7. 

2. Basic Equations for tidal perturbation 
and apsidal motion for nearly circular orbits 

Consider a binary system where one of the stars, taken to be a point object, has mass 
M and its companion star, the secondary, has mass M*. The position of the primary star 
in the orbit, as seen by the secondary star, is specified by -R(t), the separation between 
them, and the angular position, (j)o{t). The secondary is assumed to be rotating as a solid 
body with angular speed about an axis perpendicular to the orbital plane (taken to be 
the z-axis). 

The perturbation to the density of the secondary due to the tidal gravitational field 
of the primary can be written as a sum over its normal mode amplitudes: 



Sp{r,t) = J2Anim{t)Spni{r)Ye^{0,(l)), (1) 

n£m 

where Sp^i is the normalized density eigenfunction, Y£m(^,</') is the spherical harmonic 
function, and Antm is the mode amplitude which is calculated from the tidal gravity of the 
primary as described below. The pressure (p) and gravity (g) modes are uniquely specified 
by three numbers: n, the number of nodes in the radial direction, the spherical harmonic 
degree, ^, and the azimuthal number, m. We use a short hand notation, a, to denote the 
collective index (n,^, m). For conciseness of notation, we take n to be positive for p-modes 
and negative for g-modes. 
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The gravitational potential, at a point R(t), due to the tidally perturbed star is: 

^(R,t) = -G I dr ^^""^ 



|R-r|' ' ' 

Substituting the normal mode expansion of the density from equation (1) and using the 
standard spherical harmonic expansion of 1/(|R — r|) (Jackson, 1962), we obtain the fol- 
lowing expression for the perturbation to the gravitational potential: 

,T, (T> 4\ A ^ AnlmYlmjOo = Tr/2, (f>o{t) - flj)Qnl , . 

^i(R,t) = -4nG^ (2/ + ' 

nlm 

where the orbital plane is taken to be 6o=tt/2^ 

Qni= / drr'+^Sp^ir) (4) 
Jo 

is the overlap integral for the mode, i?* is the radius of the star, and the displacement 
eigenfunctions, are normalized to have unit energy, i.e.. 



(5) 



The equation for the mode amplitude, Aq,, can be cast in the form of the following 
second order ordinary differential equation for a forced, damped harmonic oscillator: 

^ + 2T.^+ulA. = Ut)^ (6) 

with 

Ut) = (2^+1) Y.^i^o = 7./2,Mt) - ^.t) = f,^ (2^^1)a^+i' 

as derived for a nonrotating star by Press and Teukolsky (1977), where Fq, is the dissipation 
rate of the mode, a is the semimajor axis of the orbit, and where the reduced forcing 
function, fim{t) = Yf*^[TT /2^ (j)o{t) — 0*t)/ [i?(t)/a] ^''""^ , contains all the time dependence 
of fa{t). The amplitude is as seen from a frame rotating with the star. The rotating 
coordinate system is related to the non-rotating coordinate system through a rotation 
about the z-axis by an angle of 0*t. 
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The solution of equation (6), in terms of the Green's function, is 
A^t) = A4t)j-—/—— = ^ / dhexp[-T4t-h)]sm[u'^{t-h)]f,^{t), (8) 



where to'^ = ycu^ — to a and where Aa{t), the reduced mode amphtude, contains aU 

the time dependence of Aa{t). 

We next derive an expression for and the periastron advance for nearly circular 
orbits, and then present results for orbits of arbitrary eccentricity. For nearly circular 
orbits (j)o{t) ~ Oot, where Oq is the angular velocity of the star, and R{t) is approximately 
constant. With these two simplifications the integral in equation (8) is straightforward to 
carry out and the resultant mode amplitude is given by: 

Lo'^ — m^O^ — i2mTa^i ' 

where 

Oi=Oo-0* (10) 

Inserting this into equation (3) and performing the sum over m=0 and m = ±2 for the 
dominant £=2 modes, we find for the perturbation to the gravitational potential. 



-ttG^Mv-^. r 3(1-402/^2 



5i?6(t) 



12 



1 + 



a / 



40?M)2 + i6F202/^4 



'11^ 



We have ignored the splitting of the eigenfrequencies due to rotation in the above equation. 
The perturbation to the gravitational potential of a star due to tides, in the static tidal 
approximation, is given by (Schwarzchild, 1958): 

-2kRlGM 

= ' (12) 

where A;, the dimensionless apsidal motion constant, depends on the structure of the star. 
Comparison of equations (11) and (12) leads to the following expression for the apsidal 
constant k. 



3(1-40?K: 



40?M)2 + 16F2 02/^4 



'131 
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In the limit that |0i | <^ Ua for the low order modes, which have the largest values of the 
overlap integral, the star can adjust its shape to the instantaneous tidal force, i.e., the tide 
is static, and the above expression for k reduces to 

* n 

The classical apsidal formula uses the value of k given by the above expression. Thus, we 
expect deviation from the classical formula (see eq. [13]) when the frequency of one of the 
low order modes is comparable to 20i. This result, as we show in the next section, applies 
to all orbits (not just nearly circular ones), although for eccentric orbits resonances with 
harmonics of the orbital period must also be considered. We also note that the periastron 
advance is linearly related to the apsidal motion constant and thus our correction to the 
apsidal motion constant is also a linear correction to the periastron advance. 

3. The effect of resonance on the apsidal motion constant 

In this section we consider the effect of orbital resonance on a given mode's contri- 
bution to the periastron advance. In the next section (§4), we consider whether or not 
this resonance has an effect on the apsidal constant for a polytropic star by considering 
its complete normal mode spectrum. The two parameters defined below are of central 
importance in determining when the classical apsidal motion formula is valid: 



where a is the semi-major axis, Mt is the total mass of the system, 0^, the angular velocity 
at periastron, is given by: 



'161 



and e is the orbital eccentricity. For a nearly circular orbit, the two parameters, and 
v Q/ , are approximately the same, but for eccentric orbits they are quite different. As we 
show below, the apsidal motion can be quite different from that given by the static tide 
approximation when ba is small (6q, < 7). 

To start our analysis, we separate out the dependent term in the reduced forcing 
function, /^m(^,0*) = exp(zmO*t)/^m(^, O* = 0). Without rotation, the forcing function 
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is periodic with period T and we can decompose it in terms of its Fourier coefficients. 
Therefore, 



fem{t,^*) = exp(zmO*t) 



Ci^"^ sin(nOot) + ^n"^ cos(nOot) 



.n=l 



n=0 



(17) 



It is straightforward to show that Im^Dn"^^) = Cn"''= for m=0 modes. In addition, 
for n > Op/Oo, Re{Di^"^^) — sign(m)/m((7i'^'"'') for m = ±2 modes. The m = 
Fourier coefficients decrease monotonicaUy with n, whereas the m = ±2 coefficients have a 
maximum at n 20p/Oo and are larger in magnitude than the m = coefficients for n > 
Op/Oo- For n greater than about 50p/Oo, the D^^'^^ fall off as exp( — 1.3ne~°'^^Oo/f^p)/ (1 — 
e)*^^/^-* and the D^^^^ are larger than the D^^'^^ by about a factor of exp(2.6e~'^^). 

Inserting the Fourier series expansion of fim (eq. [17]) into equation (8) and carrying 
out the integral, the amplitude for mode a, in the rotating reference frame, can be written 
in terms of the above Fourier coefficients as follows: 



exp(zmO*t) 
" 2 

oo 
n=l 



E 



[Dn"^^ + iCn"^^) exp{-inflot - i(t>n) 



[D 



(tm) 



n 



ms 



l212 



n 



ms 



+ 



2D 



(tm) 



m?s'^ + 2imsdc 



where 



dr 



tan c/)^ 



2da(ms 



n 



tan </)'^ 



2da{ms + n) 
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— (^n — ms )^ " — (n + ms)^ 

It is clear from equation (18) that m=0 mode amplitudes are unaffected by rotation 
(provided that we ignore the perturbation to the mode's eigenfuntion) and are resonantly 
excited when Tq, is an integer. For nonzero m, resonance occurs when Tq, = |n it ms\. The 
magnitude of the resonant term in the series, ignoring for the moment the effect of rotation, 
is approximately C^^^J^^ Sr~^ , where [r^J is the nearest integer to Tq, and Sr-Q. = Tq, — [r^J. 
Similarly, the magnitude of the n-th non-resonant term in the series is approximately 
Cn"^^ /r-Q.- Thus, the contribution of the resonant term to the mode amplitude is of the 
same order as the sum of the non-resonant terms provided that 



E — 

n=l 



exp[ 



■1.36. 



.25^ 



= )(3/2)^, 



(20) 



Using the expression for Cn"^^ given above, we find that the resonant term is of order 
the sum of the non-resonant terms provided that ba ^ — e"*""'^^ ln(^rQ,). The evolution of 
the orbit limits the value of Sva to be greater than some minimum value, Sr^im which is 
of order 10~^ for close binaries; Sva ^ Srmin is unphysical because changes in the orbital 
period, as a result of orbital evolution, cause the phase of the mode at periastron to drift by 
TT in Sr^\^ orbits (see Kumar, Ao, & Quataert, 1995). Thus, resonances have a significant 
effect on the mode amplitude, the perturbation to the gravitational potential, and the 
apsidal constant only when ba is small (less than about 7 for Srmin ~ 10~^). 

For our subsequent discussion of the validity of the classical apsidal motion formula, 
it is useful to define the following dimensionless number, which gives the fractional 
difference between the gravitational potential, due to a mode of order n, in the static tide 
limit and in the exact calculation: 

Z-^Z 1,1 

where is the perturbation to the gravitational potential arising from the n-th mode 

and ^ is the perturbation to the gravitational potential in the static tide limit (obtained 
by setting Aa = fal^lt in eq. [3]). Using equations (3) and (8) to determine and 
equations (12) and (14) to determine we obtain: 



Cn 



Qn2 



2 



^''^^l"''" V i.2m(t = 0)r2m(^o = 7r/2,</>o = 0) 

m=0,±2 



(22) 



where t = corresponds to the periastron and Rperi = a(l — e) is the periastron distance. 
By definition, C„ = in the static tide limit, whereas (n = 0.5, for example, would mean 
a 50% deviation from the classical formula because of a resonance involving a quadropole 
mode of order n. 

Since the tidal potential falls off as 1/ , the periastron advance is primarily deter- 
mined by the tidal potential at periastron. That is why we have chosen the mode amplitude 
at t = (periastron) in equation (22). We note, however, that (n is only an approximate 
characterization of the deviation from the classical apsidal motion formula. The correct 
and accurate description of the deviation from the classical apsidal formula is provided by 
explicitly calculating the periastron advance using the perturbation to the gravitational 



potential, (eq. [3]), and comparing the result with the classical formula. We thus in- 
troduce a parameter rjm which gives the fractional deviation, arising from the n-th mode, 
from the total periastron advance given by the classical apsidal motion formula. 

We find from numerical calculations that the functional dependence of (n on the modal 
and orbital parameters is very similar to that of rjn- Figure (la) shows the comparison 
between (n and function of ba- The stronger ba dependence of rjn is due to 

the stronger ba dependence of the mode amplitude away from periastron (Kumar, Ao, 
& Quataert, 1995). The discrepancy between rjn and (n is most significant for ba ^ 1. 
Since we are primarily concerned with deviation from the classical apsidal formula in the 
neighborhood of a resonance, figure (lb) shows the functional dependance of (n and rjn 
on Sva- Because (n and rjn are roughly equal, as illustrated in figure (1), we present our 
subsequent analysis of the validity of the classical apsidal motion formula in terms of the 
more analytically tractable quantity, (n- 

In this section, as we are only considering individual modes and their contribution to 
the apsidal motion formula, we take QjS = Sa^/sQa or Q\l^jjQ'^jj = 1 and investigate in 
some detail the effect of a resonance on C„ as a function of orbital parameters. This drastic 
simplification of the normal mode spectrum to a single mode is made to gain physical 
insight. In the next section we carry out the consistent calculation and include the entire 
modal spectrum of polytropic stars. We now substitute the Fourier series expansion of the 
mode amplitude (eq. [18]) into equation (22). We note, however, that the Fourier series 
technique employed in this paper is only applicable so long as the orbital evolution time 
scale is such that the resulting change in Va per orbit is less than {SraY (this condition is 
obtained by requiring that the tidal forcing of a mode is in phase at periastron over IjSra 
orbits, in order for the the mode amplitude to build up). Keeping only the d = 2 terms in 
equation (22), we obtain 







where 






and 



ni = n + 2s, n2 = n — 2s, n3 = n. 



(25) 
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The sum over j in equation (23) corresponds to the sum over m = 0, ±2. Note that we 
have again ignored the rotational splitting of the eigenfrequencies in the above derivation. 



For most physical cases of interest, we can accurately take Re{Dn ) ~ —Im{Cn ) 
and thus Ai 0. We see from equation (23) that resonance occurs when Tq, = |n — 2s| 
or when Tq, = n, causing the mode's contribution to the apsidal motion to differ from its 
static tidal contribution. 

The resonant term in equation (23) can be written as {0.5Sra /va) / [Sr'^ + c/^). Thus, 
for da ^ S^a, Cn scales as <i~^. For da ^ ^ra, C„ is independent of the damping, 
and is proportional to 5r~^ . In addition, for < 0, (n can be less than zero, which 
corresponds to apsidal recession, a result opposite to that expected from the static tidal 
formula. The dependence of (n on Sva is shown in figure (2) for several values of 6q,; 
is taken to be zero and da = 0.0005 < Sva- It is clear that (n scales as Sr~^ for small 
Sva- However, for larger Sva^ when the off-resonance terms in the Fourier series make an 
increasingly dominant contribution to the periastron advance, (n falls off more slowly. For 
&a ~ 8, C„ does not depart much from zero because of the small value of the resonant 
term that is limited by a finite damping ((/„) as well as the lower limit to |^rQ,| imposed 
by orbital evolution. 

The effect of a resonance associated with m = ±2 modes on (n depends on the stellar 
rotation rate, whereas m = resonances are, to first order, unaffected by stellar rotation. 
From equations (23)-(25) we see that, for ba ^ 20* /Op, the effect of stellar rotation is to 
increase the resonance frequency by 2s, and thus the effective ba value of the mode, near 
resonance, by about 2sOo/Op = 20* /Op. 

For ba ^ 20* /Op the dependence of (n on the stellar rotation rate is more complicated 
due to an additional possible resonance at the nearest integer to 2s — Tq,, for which the 
effective ba value of the mode is (20* — Ua)l^p- Thus if O* < Ua < 20*, rotation causes 
the effective ba value of the mode to decrease further below 2 causing to decrease. 
However, for Ua ^ O*, the effect of rotation is to increase 

For most physical systems, it is highly unlikely to find a low order p ov g mode 
with ba ^ 20* /Op and thus the effect of rotation is to increase Ua by about 20* and 
consequently to decrease However, the quasi-toroidal modes that exist in rotating 

stars have Ua ~ O*, and thus their effective ba values, near resonance, can be about 2, 
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thus enhancing their contribution to the periastron advance (this is discussed in greater 
details in §5). 

Finally, we note that for a fixed 6q,, (n has a relatively weak dependence on the orbital 
eccentricity, for e > 0.2. In figure (3) we show (n as a function of ba for e = 0.4, 0.6, and 
0.8. Note that different values of the eccentricity, for a fixed value of 6q,, correspond to 
different orbits with the same periastron passage time. Physically, the weak eccentricity 
dependence is due to the fact that the perturbation to the gravitational potential due 
to tides falls off rapidly away from periastron. Thus, for eccentric orbits, most of the 
contribution to the apsidal motion comes from near periastron and so orbits of different 
eccentricity but the same periastron passage time have nearly the same (n- The more 
rapid decrease of (n with ba for smaller e arises from the strong eccentricity dependence 
of the Fourier coefficients for small e (recall that oc exp( — 1.36Q,e~°'^^)). For e < 0.2, 
the dominant resonance is at 2, and this case is described well by equation (13) 

of the previous section. 



Having discussed in some detail the effect of a resonance on the contribution of an 
individual mode to the apsidal motion, we now apply this discussion to the realistic sce- 
nario of a star with some prescribed normal mode spectrum. In particular, we consider 
polytropic stars and investigate when to expect significant departures from the classical 
apsidal formula. 

It is useful to express the mode frequencies of polytropic stars in the following dimen- 
sionless form 



This nondimensional frequency is a function of polytropic index (upi) but is independent 
of stellar mass. We rewite the parameter 6q,, defined in equation (15), as follows: 



where dmin = a[l — e)/i?* is the distance between the two stars at periastron in units of 
the stellar radius. 



4. Apsidal Motion for polytropic stars 






.1/2 



(27) 
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As was shown in the previous section, the classical apsidal motion formula can be in 
error when the value of ba for at least one of the low order modes of the star is less than 
some small value, 6c, which is of order 7. This condition, together with equation (27), 
implies that the classical apsidal formula is applicable so long as 

d^^n>[j^j [ \ ■ (28) 

It should be emphasized that ba < be is a necessary, but not sufficient, condition for the 
breakdown of the classical formula. Consequently, equation (28) gives a sufficient condition 
for the validity of the classical apsidal formula. 

The /, gi and pi modes of stars generally have the largest overlap integrals (IQnH), 
and since k is proportional to Q^i, they contribute the most to the apsidal constant. 
Therefore, if a resonance associated with a mode is to cause deviation from the classical 
apsidal formula, its contribution to the apsidal constant must be greater than approxi- 
mately {Qo2 / Qn2)'^ times its classically expected value. 

The overlap integrals decrease rapidly with the number of radial nodes, n, as is shown 
in figure (4) where we plot the overlap integrals, in units of (i?^/G)"^/^, as a function of 
n for several polytropes. The decrease of |(5n2|, for fixed n, with Upi does not imply that 
stars of higher polytropic index are more likely to lead to deviation from the classical 
apsidal formula since mode frequencies and thus ba values also increase with Upi. The 
rapid decrease of \Qn2 \ with n means, however, that deviation from the classical formula is 
most likely to arise due to a resonance involving the / or low order g modes. Resonances 
involving higher order g modes cannot modify the classical result much, unless the modes' 
overlap integrals happen to be within a factor of a few of the /-mode. 

We now quantify the likelihood of significant error in the classical apsidal motion 
formula due to resonance involving low order modes. Since the orbital and modal frequen- 
cies are independant quantities, Sva should be treated as a random number with a value 
between -0.5 and 0.5. Thus, the value of |^rQ,| required to get a certain value of (n is a 
good measure of the probability of failure of the classical apsidal formula. We thus define 
a probability function, which is equal to twice the largest value of |^rQ,| such that a 
resonance involving a mode of order n leads to a fractional error of (n- Pn is a function of 
6q,, c/q,, s, and the stellar structure. The dependence of on ba and Qni is determined 
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by separating out the resonant term in equation (23) and solving for Sva in the limit that 



P„ OC .291 



CnZ^Ql/Ql + ^l -9{ba) 

where Dr^ is the Fourier coefficient of the resonant term in the Fourier series and g{ba)^ 
the sum of the non-resonant contributions to is a generally increasing function of ba- 
Thus, for modes that make a small contribution to the periastron advance in the static 
tide limit, i.e., Qj/Q^ ^ 1, or for large we have a simple scaling law for 

In figure (5), we show the results of calculating numerically for several different modes 
and polytropic indices. For all of these calculations we have taken Tq, = 15.0, da = 0.0005, 
and = 0. Note that the main features seen in figure (5) are just what we expect from 
equation (30). For ^ 0, the effect is simply to shift the ba values, as discussed in §3. 
This is explicitly illustrated in figure (6) which shows as a function of ba for the / 
mode of an Upi = 3 polytrope for three different rotation rates. In addition, the weak Tq, 
dependence of P^ follows from the dependence of the Fourier coefficients as discussed in 
the previous section. This is because calculations of P^ for two different values of Tq, and 
eccentricity, for a fixed ba, are related by the ratio of the Fourier coefficients corresponding 
to their resonant terms (see eq. [30]). Finally, provided that da ^ Sr^in (as is likely to be 
the case), Sr^im not da, is the smallest physical value of Sva, and thus P^ is independent 
of the dissipation rate. 

The probability of deviation from the classical formula, P^, due to a resonance in- 
volving higher order g and p modes of these polytropes can be obtained from figure (5) by 
scaling the variables according to equation (30) and using the overlap integrals from figure 
(4). For example, the probabilities of obtaining (n = 0.3 due to a resonant gs mode of an 
Upi = 2.0 polytrope can be calculated using equation (30) and figure (4), and we find them 
to be about 20 times smaller than those corresponding to (n = 2.0 for the gi mode of an 
Upi = 3.0 polytrope which are shown in part (c) of figure (5). 

Finally, we illustrate how to use the above results to explicitly estimate the system 
parameters for which the apsidal motion formula should be applicable. For example, to 
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be 90% certain that the static tidal result is within 50% of the exact result, i.e. (n < 0.5, 
for a npi = l.b polytropic star, requires that the ba value of the / mode should be greater 

1 /3 

than about 5 (see fig. 5d). This means that for dmin > 2.2 [(1 + e)MT/M^] ' , we expect 
the classical apsidal formula to be valid to the degree specified. 

5. Effects of rotation modes on apsidal motion 

In a rotating star, the presence of the coriolis force provides an additional restoring 
mechanism for oscillatory motions. Thus there exists a new family of modes, called the ro- 
tation modes or r-modes, which have frequencies of order the stellar rotation rate. Because 
of their small frequencies, r-modes are likely to be resonantly excited by tidal interactions. 
In this section we examine the effects of r-modes on apsidal motion. 

In a non-rotating star, purely toroidal modes with V • ^=0 and (^r=0 exist as part 
of the general solution for nonradial oscillations. These modes correspond to horizontal 
displacements of the fluid associated with which there is no pressure perturbation; there- 
fore, their mode frequencies are zero. Moreover, these modes are not excited by the tidal 
force since their density perturbations are zero. However, the toroidal modes become non- 
trivial in a rotating star. Papaloizou and Pringle (1981) and Rocca (1982) have examined 
in some detail the importance of r-mode resonances in the context of tidal excitation. We 
closely follow the approach of Provost, Berthomieu, and Rocca (1981) (hereafter, PBR) to 
calculate the eigenfunctions and frequencies of r-modes in a slowly and uniformly rotating 
star. We carry out a perturbative expansion to the structure and eigenmode equations 
of the star using the non-dimensionalized rotation rate = Q^R^/ GM^ as a small ex- 
pansion parameter. We also assume, as did Papaloizou and Pringle and PBR, that the 
perturbations are adiabatic and that the perturbation to the gravitational potential can 
be neglected (Cowling approximation). Saio (1982) has shown that both of these approx- 
imations are quite reasonable. Our primary interest lies in the coupling of the r-modes 
to the quadrupole gravitational potential, and for this we need to calculate their overlap 
integrals. The basic equations are given in the Appendix. In the following, we summarize 
the salient features of the r-modes, as relevant for our purposes. To lowest non- vanishing 
order in O*, the eigenfrequencies are given by: 
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2m 

'^ntm = J^^p^'^^*^ (31) 

For m = 0, cUq,=0, and thus there are no axisymmetric r-modes. To this order, the 
r-mode eigenfunctions correspond precisely to the toroidal modes of a nonrotating star 
described above; therefore, even modes with m ^ Q are not excited by the tidal force. 
Hence we must go to the next order. Higher order physical quantities, as it turns out, 
can be expressed in terms of lower order quantities. The spherical harmonic functions, 
however, can no longer adequately describe the rotationally distorted star. The next 
order quantities have an angular dependence described by a superposition of two spherical 
harmonics: Y^+i^m, and Yf-i^^- This means that one has to be careful in calculating the 
overlap integrals from equation (4) since in its derivation the angular integration has been 
performed under the assumption that the mode eigenfunctions are well described by Yf^- 
Since the eigenfunctions now consist of sums of Yi^i^m and Y^-i^^, it follows from the 
orthogonality condition of the spherical harmonics that in order to calculate the overlap 
integrals (5^2, ±2, where the subscripts (2, ±2) now refer to the = 2, m = ib2)th term in 
the expansion of the gravitational potential, we need the ^ = 3, m = ±2 r-modes (instead 
of ^ = 2, m = ±2 for the p and g-modes, which are based on a spherically symmetric star). 
We also find that to lowest nonvanishing order, Qa oc O*. It is also worth noting that the 
higher order correction to the eigenfrequency varies as which, unlike the lowest order 
term, depends on the order of the mode. We refer the reader to the Appendix for more 
details. 

In previous sections, we have shown that the tidal excitation of a mode is primarily 
determined by two parameters, ba and Qa- We first discuss the ba dependence of the 
r-modes. Equation (31) shows that for ^ = 3,m = ±2, the lowest order frequency is given 
by 0*/3. Note that this is independent of the number of radial nodes (n) because, to this 
order, the mode is purely toroidal. Ignoring higher order corrections, we find ba = 0*/30p. 
However, as shown in §3, stellar rotation effectively shifts the ba value of these modes 
(near resonance) by 2(0* —LOa)/^p so that b^f^ = (5/3)0*/Op. For synchronous rotation 
at periastron, b^f^ = 5/3, which is quite favorable for departure from the classical apsidal 
formula since m = ±2 modes couple best to the tidal force for b^f^ 2. Non-synchronous 
rotation moves b^f^ away from its optimal value. Thus, (n decreases exponentially for 
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super-synchronous rotation. 

Using equation (A6) in the Appendix, we compute the quadrupole overlap integrals 
for several polytropes [upi = 1,2, and 3). The numerical results for the low-order modes 
are shown in Table (1). Note that the values in the table are given in terms of the 
nondimensionalized overlap integral Qa (Qtm = Qa sj R%l G)- Our results show that the r- 
modes are unlikely to make a significant contribution to the apsidal motion in spite of their 
optimal hct values. The reason for this is that their overlap integrals are too small compared 
to that of the f-mode. /,From Table (1), we see that the overlap integrals are generally 
largest for stars of smaller polytropic index. This is to be expected since stars of higher 
polytropic index corresponds to more centrally condensed stars, and the effects of rotation 
for these stars are generally smaller. If, for example, =0.1 for a Upi = 1 polytrope, the 
overlap integral for the n = r-mode is three orders of magnitude smaller than that of 
the f-mode. This makes it unlikely to cause significant deviation from the classical apsidal 
formula. The situation is not much improved for a more rapidly rotating star since the 
overlap integral increases only linearly with O*. Indeed, even if we extrapolate our results 
to the break-up rate (i.e. = 1), the overlap integral for the lowest order r-mode of an 
Upi = 1 polyrope is still about one order of magnitude smaller than that of the f-mode. 
Nevertheless, near break-up rotation, we are no longer justified in ignoring higher order 
corrections, and our results could be significantly modified. 

To summarize, we find that the r-modes are unimportant for the apsidal motion of a 
slowly and uniformly rotating star. The same conclusion is expected to hold for rapid, but 
non-synchronous, rotation. This is because the overlap integrals are expected to grow only 
as some power of the rotation rate (even when we include higher order corrections) while C„ 
decreases exponentially with increasing ha- Hence, the only regime where r-modes might 
be important is for near break-up, synchronous rotation, but other methods are needed to 
investigate this regime in detail. 

6. Application to AS Cam and DI Her 

Close binary systems offer a unique opporunity to test stellar structure theory as well 
as provide tests for various theories of gravity. Two binary systems, AS Camelopardalis 
and DI Herculis, are well known because of their anomolous apsidal motions (Guinan and 
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Maloney, 1985; Maloney et al. 1989; KhaliuUin et al. 1983, 1991) which are, respectively, 
about 1/3 and 1/7 of the classically expected values. In this section, we investigate whether 
resonances of stellar modes with the orbit might be responsible for their anomolous be- 
havior. 

AS Cam consists of two main sequence stars (Mi = 3.3Mq^ Ri = 2.^7 Rq, M2 = 
2.5M0,i?2 = I.9-R0) in an e = 0.1695 orbit with a period of approximately 3.4 days. 
The damping time of low order modes of main sequence stars is of order 10^ to 10"^ years; 
therefore, da ~ 10~^ . Main sequence stars in this mass range have typically both convective 
and radiative zones and thus are not accurately represented by a single polytropic index. 
As a crude aproximation, however, we take the stars to be modelled as Upi = 2.5 polytropes 
and calculate the mode frequencies and overlap integrals. 

In previous sections, we have, for simplicity, assumed that one of the stars was a point 
mass. This restriction can easily be relaxed by applying our method to each star and 
ignoring the second order effects. Using an Upi = 2.5 polytrope, we determine that the 
value for the f mode of the first star (Mi = 3.3M0,i?i = 2.^7Rq) is about 20, whereas 
the hct for the f-mode of the second star is 27.4. Because ^ — In(^rmm), it is impossible 
that an f-mode resonance can give rise to the observed apsidal motion. A more favorable 
value of ba is obtained for a high order g-mode. For instance, the (/15 mode of the star 
has ba ~ 2.2 and / ~ 10^. In order to obtain (n = —0.66, as required to account 

for the observed apsidal motion, the n = 15 mode of an Upi = 2.5 polytrope must have 
\Sra \ ~ 10~^. This degree of resonance is, aside from being improbable, unphysical, since 
it is impossible for the mode to have phase coherence over the required 10^ orbits, that 
is, Sva ^ Srmin- Because of their higher frequencies, p-modes are even less likely to be 
resonantly excited. Thus, we conclude that the anamolous apsidal motion of AS Cam is 
not the result of resonances involving g, f or p-modes. 

The stellar rotation rates for both stars in AS Cam are observed to be 1.05Oo ~ 

0.73Op (Maloney et al. 1989). Thus the stars are rotating somewhat sub-synchronously. 
The r-modes therefore have b^f^ 1.2. From Table (1), (interpolating between Upi = 2 
and Upi = 3), we find that the nondimensionalized overlap integral for the n = mode 
(which generally has the largest overlap integral) is 0.04O*. (Recall from the previous 
section that, to lowest order, Qa oc O*.) With the quoted value for O*, 0.08 (for star 
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1), implying that - I Q\ ~ 10^. Thus, scaling by equation (30), we find |^rQ,| ~ 
to account for the observed apsidal motion via an r-mode resonance. We can therefore also 
exclude r-mode resonances as a possible explanation for the anamolous apsidal motion of 
AS Cam. 

DI Herculis consists of a primary with Mi = S.lM© and Ri = 2.68Rq and a secondary 
with M2 = 4.5M0 and R2 = 2ASRq in an e = 0.489 orbit with a period of about 10.55 
days. We also use an Upi = 2.5 polytrope as a crude model for these stars. Furthermore, 
since the periastron passage time and M^/R^ are larger for the DI Herculis system than 
for AS Cam, the ba values of the modes are slightly larger than those of AS Cam. Since the 
dominant dependence of is on 6q,, we again conclude that f, p, and g-mode resonances 
cannot cause the observed deviations in the apsidal motion. In addition, the stellar rotation 
rates of both stars are observed to be 3.50o ~ (Guinan and Maloney 1985). Thus 

the two stars are rotating synchronously at periastron, implying b^f^ = 1.7 for the r- 
modes. Nevertheless, for the primary star, 0.02, so that the overlap integrals are 

even smaller than those of AS Cam. Thus, as before, we can safely rule out r-modes as a 
possible explanation for the anamolous apsidal motion. 

We conclude that the anamolous apsidal motion of DI Herculis and AS Cam is not 
the result of the breakdown of the static tidal approximation. Other possible explanations, 
such as the presence of a third body or an inclined spin axis with respect to the orbital 
plane, are discussed in detail in the above references and references contained therein. 

6. Summary and Discussion 

We have investigated the validity of the classical apsidal formula for tidal distortion by 
calculating the perturbation to the gravitational potential using a Fourier series expansion 
of the mode amplitude; this expansion is a good approximation so long as the orbital 
evolution time scale is such that the change in Tq, in one orbit is less than [SraY (where 
Tq, is the ratio of the mode frequency to the orbital frequency and Sva = — where n 
is the nearest integer to Tq,). 

We find that the classical apsidal motion formula, which assumes static tides, is very 
accurate when the dimensionless number 6q,, which is the ratio of the periastron passage 
time to the mode period, for the low order f, p, and g-modes, is greater than about 7. 
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However, when 6q- ^ 7 for one of these modes, dynamical effects become important and 
the periastron advance can differ significantly from that computed using the static tidal 
formula. The deviation is largest when one of the low order modes is nearly resonant with 
the orbit, i.e., when the mode period, for a non-rotating star, is some integral fraction of 
the orbital period. For a near resonant mode of small ba the deviation from static tidal 
formula scales as Sr~^ , so long as the mode damping time divided by the orbital period is 
greater than Sr~^ . 

The effective ba value, near resonance, of a mode in a star rotating synchronously at 
periastron is generally larger than in a non-rotating star by about two; therefore, the static 
tidal approximation is better for the p, f, and g-modes of a rotating star. 

The contribution of a given mode to the apsidal constant is proportional to Q^i, where 
Qnh the overlap integral, is usually largest for the / mode and decreases strongly with the 
number of radial nodes. Thus, deviation from the classical apsidal motion formula is only 
likely to arise because of a resonant / mode or a resonant low order g mode whose overlap 
integral is within a factor of a few of the / mode overlap integral. We have quantified 
the probability of deviation from the classical apsidal formula by calculating how close the 
frequency of a low order mode must be to a resonance, Sva^ in order to obtain a certain 
deviation from the classical formula as a function of ba for several polytropes (§4). 

For rotating stars there exists a new family of modes (r-modes), where the restoring 
force is the coriolis force, which have frequencies of order the stellar rotation rate. These 
modes have optimal effective ba values and can easily be resonantly excited by the tidal 
force. We have investigated their effect on apsidal motion for slowly and uniformly rotating 
stars and find that their overlap integrals are generally too small for them to be important. 
Only when the star is simultaneously rotating near break-up and synchronous at periastron 
could the r-modes possibly become important. 

We have applied our dynamical calculations of the periastron advance to AS Cam and 
DI Herculis in order to determine if the discrepancy between the observed and predicted 
apsidal motion for these systems is due to the break down of static tidal approximation 
for these systems. We conclude that this is very unlikely to be the case. Given the 
observational failure to find a third body in either system, we suspect that the anomalous 
apsidal motion for these systems is likely to be caused by an inclined spin-orbit coupling. 
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as has been suggested by Shakura (1985) and Company et al. (1988). 
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Appendix 



In this Appendix, we summarize the basic equations for r-modes in a slowly and 
uniformly rotating star as described in PBR. We also give the expression for the overlap 
integral. 

Express all physical quantities in the natural units where = i?* = G = 1. Let 
be the rotation rate of the star and let the z-axis be the axis of rotation. We work in the 
rotating frame of the star. Since is assumed to be small, we can express the equilibrium 
structure of the st 

r = x{l -nla{x)cos^ e), (Al) 

where x labels the equipotential surfaces, and a[x) describes the rotational distortion to 
lowest order. We calculate a[x) numerically using the Clairaut-Legendre expansion as dis- 
cussed in Tassoul (1978). Similarly, we can expand the eigenfrequencies and eigenfunctions 
in powers of as follows: 

(jj = O*c<;o(l + O^cui) 



The zeroth order fluid equations correspond to the conservation of mass and angular 
momentum about the z-axis. Its solutions are 



2m 







imCimi 



Pr(cos^) 



sin 6 



■,-C,^(x)^Pr(cos^) 



(A3) 



where Cim{x) describes the mode structure in the radial direction (unrelated to the Fourier 
coeflicients deflned earlier!). To zeroth order, Cim{x) is completely arbitrary (which ex- 
plains why the zeroth order frequency is independent of the order of the mode); it is 
determined only through the flrst order equations. PBR derived the following second 
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order ODE for Cim{x): 



— Xiuji h A— In — 



dx"^ dx \ Ag J dx i x"^ x"^ dx \ g 

Ctm = 



(A4) 



where p^P^g, and 7 have their customary meanings, and 

^ _ 1 c/p I dP 

p dx 7P dx 

The As depend only on d and |m| and are defined as: 
Ai = {).2bDit^{l+lf 

A2 = Dt[{l + lfK.{l - l)A + (^) - eK.{^)K+{l + 1)] 
A3 = {).bDif{l+lfm^ 

A4 = Q.hDd^{l+lf[{l + lfK-{l-l)K+{(l)+eK-{(l)K+{l+l)] 



wi 



ith 



2^+1 
l-\m\ + l 



21+1 

Di = [{£ + l)^K-{£ - l)K+{£) + £^K-{£)K+{£ + 1)]-^ 

It is clear from equation (A4) that Cim is independent of the sign of m. Having 
obtained Cirm one can then go back to the first order fiuid equations and solve for the 
first order quantities in terms of Cim- Consequently, the overlap integral, which can be 
expressed in terms of Cim and its derivative, depends only on |m|. Hence, for notational 
simplicity, we assume that m > in subsequent discussions. Since Cim gives the horizontal 
displacement, it must be zero at the center and finite at the surface. Equation (A4), 
together with these boundary conditions, constitute a well-defined eigenvalue problem 
with eigenvalue ui. We have numerically computed C32{x) and the corresponding ui for 
several polytropes since these are the modes that couple to the tidal force (see below). 
The results for toi are shown in Table (2). Note that when the correction frequency toiQl 
becomes comparable to 1, equation (A4) is no longer valid. Since |c<;i| tends to increase 
with mode order n, for a given O*, there are only a limited number of low to moderate 
order modes we can obtain from the PER procedure. 
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To calculate the overlap integrals, we start with a definition more general than equa- 
tion (4): 



Qn 



£m 



£'m' 



(A5) 



where the factor of Yf*^ comes from the standard expansion of the gravitation potential 
1/|R — r| in terms of spherical harmonics. It is easy to see that the above definition reduces 
to equation (4) when the displacement eigenfunction involves only one spherical harmonic. 
From the above equation, we can readily calculate the overlap integrals in terms of Cim 
and obtain, to lowest non- vanishing order. 
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We normalize the eigenfunctions such that each mode carries unit energy. Hence 
Cim 1/0* and Qa oc O*. Note that in general, to calculate the overlap integral cor- 
responding to the (^, m)th term in the expansion of the gravitational potential, we need 
the r-mode eigenfunctions labelled by [£ -\- l,m) and [£ — l,m). For tidal excitation, we 
are mainly interested in quadrupole modes with ^ = 2, m = 2, in which case we need to 
calculate C32{x) to obtain the required overlap integrals. 
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Figure Captions 



FIG. 1. — A comparison of rjm the fractional deviation, for a given mode, from the perias- 
tron advance predicted by the classical apsidal formula, with the fractional deviation, 
for a given mode, from the perturbation to the gravitational potential at periastron pre- 
dicted by the classical apsidal formula. Part (a) of this figure shows (n and rjn as functions 
of ba while Part (b) compares the two as functions of Sva- Since most of the contribution 
to the periastron advance comes from the perturbation to the gravitational potential near 
periastron, as illustrated by this figure, is a reasonable approximation to rjn] conse- 
quently, in this paper, we adopt (n as an approximate measure of the deviation from the 
classical apsidal formula. 

FIG. 2. — The fractional deviation (Cn), for a given mode, from the classical apsidal formula 
as a function of Sva for different values of ba- = 0.0005 < and thus C„ oc 5r~^ . 
This dependence, however, sets in at increasingly larger values of 8rct for increasing values 
of bet as the non-resonant terms in C„ make an increasingly dominant contribution to the 
periastron advance. By b^ ~ 8, the 5r~^ dependence only sets in at ~ and thus the 
resonance has little effect. We therefore emphasize that resonances effect C„ and deviation 
from the classical apsidal formula is only expected for small b^- 

FIG. 3. — C,n as a function of b^ for three values of the orbital eccentricity: e = 0.4, 0.6, 
and 0.8. Keeping e fixed and varying b^ corresponds to varying r^- For a fixed value of 
6q,, the deviation from the classical apsidal formula is approximately the same for orbits 
of different eccentricities. Physically, this is due to the fact that the tidal potential falls 
off rapidly from periastron and thus the coupling of a given mode to the tidal force is 
primarily determined by the periastron ee time. 

FIG. 4. — The overlap integral, Qn2-, in units of as a funtion of the number of 

radial nodes, n, for stars of different polytropic index, Upi. n < is for g modes and n > 
is for p modes. The value of the apsidal constant is proportional to and thus the 

/ mode and low order p and g modes contribute the most to the apsidal constant. The 
rapid decrease of Qn2 with n therefore implies that high n g modes, although they are 
more likely to have 6 1, are unlikely to cause deviation from the classical apsidal formula 
because of their small overlap integrals. 
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FIG. 5. — Pm defined to be twice tlie largest value of Sva such that a certain value of (n is 
obtained as a function of 6q,, s, etc., is a good measure of the probability of failure of 
the classical apsidal formula. Parts (a), (b), and (c) of this figure show for several values 
of (n as a function of ba for the /, (/i, and pi modes of an Upi = 3.0 polytrope, respectively. 
Part (d) shows P^ for the / mode of an Upi = 1.5 polytrope. We find, however, that the / 
modes of polytropes with indices between 1.5 and 2 exhibit the same deviation from the 
classical apsidal formula so these results are also valid for Upi = 2. P^ for additional p and 
g modes of these polytropes and for additional (n can be obtained from the results in this 
figure by a simple scaling law given in equation (30). 

FIG. 6. — Pm for (n = 0.5, as a function of ba for three different rotation rates: = 0.0, 
= 0.25Op, and = O.50p, where is the angular velocity of the star at perias- 
tron. For modes with 6q, > 1, rotation effectively increases the ba value of the mode by 
approximately 20* /Op, thus making the static tidal approximation more applicable. 
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